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INFLUENCE  OF  STRONG  SELF  ELECTRIC  FIELDS  ON  THE  ION  RESONANCE  INSTABILITY 


IN  A NONNEUTRAL  PLASMA  COLUMN 

Ronald  C.  Davidson* 

Division  of  Magnetic  Fusion  Energy 
Energy  Research  and  Development  Administration 
Washington,  D.  C.  20545 

Hwan-sup  Uhm 

Department  of  Physics  and  Astronomy 
University  of  Maryland,  College  Park,  Md.  20742 

The  influence  of  strong  self-electric  fields  on  the  ion  resonance 

instability  is  examined  for  a cylindrical  nonneutral  plasma  column 

Immersed  in  a uniform  axial  magnetic  field  Bq£z*  The  analysis  is 

carried  out  within  the  framework  of  a macroscopic  cold-fluid  model, 

and  electrostatic  stability  properties  are  investigated  for  the  case 

of  rectangular  electron  and  ion  density  profiles.  The  parameter 
2 2 

6*(2u>  /u)  ) (1-f ) is  introduced  as  a convenient  measure  of  the  relative 

pe  ce 

strengths  of  the  equilibrium  self-electric  force  and  the  magnetic  force 

on  an  electron  fluid  element.  (Here,  di  is  the  electron  plasma 

pe 

frequency,  u>ce  is  the  electron  cyclotron  frequency,  f«ni/ng  is  the 
fractional  charge  neutralization,  and  <5“1  corresponds  to  the  maximum 
allowed  charge  density  for  radial  confinement  of  the  equilibrium 
configuration.)  An  important  conclusion  of  this  study  is  that  the 
equilibrium  self-electric  field  can  have  a large  influence  on  stability 
behavior.  In  particular,  stability  properties  for  6<1  differ 
substantially  from  those  obtained  when  6<<1.  Moreover,  for  a 
nonneutral  plasma  column  with  significant  charge  neutralization, 
it  is  found  that  the  fundamental  mode  (£“1)  is  not  the  most  unstable 
mode.  Rather,  higher  harmonic  perturbations  have  larger  growth  rates. 
However,  in  the  limiting  case  where  f<<l  and  6<<1,  the  t»l  mode  can 
have  the  largest  growth  rate,  which  is  consistent  with  the  result 
previously  obtained  by  Levy  et  al. 

*0n  leave  of  absence  from  the  University  of  Maryland,  College  Park,  Md. 


There  have  been  several  recent  experimental  and  theoretical 

studies  of  the  fundamental  equilibrium  and  stability  properties  of 

nonneutral  plasmas'*"  in  both  mirror^  ^ and  uniform^  ^ magnetic  field 

configurations.  Perhaps  one  of  the  most  basic  instabilities  that 

characterizes  a nonneutral  plasma  with  both  electron  and  ion 

11-13 

components  is  the  ion  resonance  instability.  In  cylindrical 

geometry  (Fig.  1),  the  ion  resonance  instability  can  be  described 
as  a two-rotating-stream  instability^  in  which  the  relative  rotation 


between  electrons  and  ions  is  produced  by  the  equilibrium  self- 
electric field  E^(r)^..  Previous  analyses'*''*' 0f  this  instability 
have  been  restricted  to  low  beam  densities  and  small  values  of 


A 

Here,  w is  the  electron  plasma  frequency,  a>ce  is  the  electron  cyclotron 
frequency,  and  f=n*j*/n^  is  the  fractional  charge  neutralization 
provided  by  the  positive  ion  background.  In  this  paper,  we  investigate 
the  influence  of  strong  self-electric  fields  on  the  ion  resonance 
instability  in  a nonneutral  plasma  column  with  rectangular  electron 
and  ion  density  profiles  (Fig.  2).  The  analysis  is  carried  out  within 
the  framework  of  a macroscopic  model,  treating  the  electrons  and  ions 
as  cold  (Tj->0)  fluids  immersed  in  a uniform  axial  magnetic  field 
Moreover,  the  stability  studies  assume  electrostatic  perturbations 
with  infinitely  long  axial  wavelength  (3/3z-0). 

It  is  useful  to  introduce  the  dimensionless  parameter 
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2u> 

6 = (l-f)  , 

a> 

ce 

which  is  a measure  of  the  characteristic  relative  strength  of  the 
equilibrium  self-electric  force  and  the  magnetic  force  on  an 
electron  fluid  element.  In  the  present  equilibrium  and  stability 
analysis,  the  parameter  6 and  the  fractional  charge  neutralization 
f are  allowed  to  span  the  range  of  values 


Q < 6 < 1 , 

and 

0 < f < 1 , 


where  6=1  corresponds  to  the  maximum  allowed  charge  density  for  radial  t 

confinement  of  the  equilibrium  configuration.  One  of  the  most  important 
conclusions  of  this  study  is  that  the  equilibrium  self-electric  field 
can  have  a large  influence  on  stability  behavior.  In  particular, 
stability  properties  for  6<1  differ  substantially  from  those  obtained 
when  6<<1  (see,  for  example.  Fig.  5).  Moreover,  for  a nonneutral 
plasma  column  with  significant  charge  neutralization,  it  is  found 
that  the  fundamental  mode  (£=1)  is  not  the  most  unstable  mode. 

Rather,  higher  harmonic  perturbations  have  larger  growth  rates  (see, 
for  example,  Figs.  8 and  9).  However,  in  the  special  limiting  case 
where  f«l  and  6<<1,  the  £,=1  mode  can  have  the  largest  growth  rate, 

12 

which  is  consistent  with  the  result  previously  obtained  by  Levy  et  al. 

[Fig.  8]. 

The  outline  of  this  paper  is  the  following.  In  Section  2,  we 
give  a brief  description  of  the  theoretical  model  (Section  2. A)  and 
summarize  the  equilibrium  properties  for  general  density  profiles 
nj(r),  j=e,i  (Section  2.B).  The  corresponding  eigenvalue  equation 

4 ay"  V ' ^ 


r~ 


that  describes  the  electrostatic  stability  properties  of  the  nonneutral 

plasma  column  is  derived  in  Section  3. A [Gq.  (13)].  In  Sections  3.B 

and  3.C,  the  eigenvalue  equation  is  solved  for  the  case  where  the 

electron  and  ion  density  profiles  are  rectangular  [Eqs.  ( 20) — ( 22) ] . 

This  leads  to  a fourth-order  algebraic  dispersion  relation  that 

determines  the  complex  eigenfrequency  u>  [Eq.  (26)].  A detailed 

numerical  analysis  of  the  dispersion  relation  is  presented  in  Section  4, 

where  stability  properties  are  Investigated  for  a broad  range  of 

plasma  parameters  “pe/wce»  etc* 

Finally,  we  emphasize  that  the  present  analysis  is  based  on  a 

macroscopic  cold-fluid  model  in  which  the  ion  (and  electron)  motion 

is  assumed  to  be  laminar.  Although  this  is  a reasonable  approximation 

when  r.  <<R'  (where  r,  . is  the  characteristic  thermal  ion  Larmor 
Li  p LI 

radius,  and  Rp  is  the  radius  of  the  plasma  column),  we  expect  signif- 
icant modifications  to  the  stability  behavior  when  rLi^Rp-  T^e 

influence  of  finite  ion  Larmor  radius  effects  on  the  ion  resonance 

14 

instability  is  currently  under  investigation  within  the  framework 
of  a hybrid  Vlasov-fluid  model,  which  treats  the  electrons  as  a 
macroscopic,  cold  fluid  and  the  ions  in  a fully  kinetic  (Vlasov) 


manner . 
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2.  THEORETICAL  MODEL  AND  EQUILIBRIUM  PROPERTIES 

2. A Theoretical  Model 

As  Illustrated  in  Fig.  1,  the  equilibrium  configuration  consists 
of  a cylindrirally  symmetric  nonneutral  plasma  column  that  is  infinite 
in  axial  extent  and  aligned  parallel  to  a uniform  applied  magnetic  field 
Bq£z.  The  mean  motion  of  the  plasma  components  is  in  the  azimuthal 
direction,  and  the  applied  magnetic  field  provides  radial  confinement 
of  the  plasma  column.  As  shown  in  Fig.  1,  the  characteristic  radius 
of  the  plasma  column  is  denoted  by  R^,  and  we  introduce  a cylindrical 
polar  coordinate  system  (r,0,z)  with  the  z axis  coinciding  with  the  axis 
of  symmetry;  r is  the  radial  distance  from  the  z axis,  and  0 is  the  polar 
angle  in  a plane  perpendicular  to  the  z axis.  The  deviation  from  equil- 
ibrium charge  neutrality  [n^(r)^n^(r) ] produces  an  equilibrium 
radial  electric  field  ^(^c)=E^(r)^r  that  influences  the  azimuthal  motion 
of  the  plasma  components. 

r 

The  following  are  the  main  assumptions  pertaining  to  the  equilibrium 
and  stability  analysis: 

(a)  Equilibrium  properties  are  independent  of  z (9/3z=0)  and 
azimuthally  symmetric  (3/90=0)  about  the  z axis.  For  example, 
the  mean  equilibrium  velocity  of  component  j (j=e,i)  can  be  expressed 

0 o a 

as  V (£)*V  (r)<^0,  where  e0  is  a unit  vector  in  the  0 direction,  and 


3vJ  /30-O-3V®  /3z. 

JO  JO 

(b)  For  present  purposes,  the  electrons  and  ions  are  treated  as 
cold  (Tj-*-0)  fluids  immersed  in  the  uniform  axial  magnetic  field  Bq^. 
Within  the  context  of  the  electrostatic  approximation  (B^B^e^  and 
^x^=0),  the  equation  of  motion  and  continuity  equation  for  each  plasma 
component  (j=e,i)  can  be  expressed  as 
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If  ”j  + 7 • <"j  V ■ 0 


where  Jg(^,t)“-V$(;|c,t)  is  the  electric  field,  rij(jc,t)  is  the  density, 
Vj(jc,t)  is  the  mean  velocity,  and  and  are  the  charge  and  mass, 
respectively,  of  a particle  of  species  j.  In  Eq.  (1),  the  spatial 
variation  in  is  neglected  (low-beta  approximation) , and  the 
electrostatic  potential  <t>(jc»t)  is  determined  self-consistently  from 
Poisson's  equation 


”4>  = 4Tje(ne~n^)  , 


where  -e  is  the  electron  charge. 

(c)  In  the  stability  analysis,  flute  perturbations  with  9/32*0 
are  considered.  All  quantities  are  expressed  as  an  axisymmetric 
equilibrium  value  plus  a perturbation,  i.e.,  ^(x,t)=ij>n(r)+6iji(x,t)  • 

'b  U *b 

Assuming  perturbations  with  azimuthal  harmonic  number  £,,  the  quantity 
6i|)(^,t)  is  then  expressed  as 


6i|)(ji,t)  = (54'j,(r)exp[i(£0-a)t)  ] , 


where  w is  the  complex  eigenfrequency. 


2.B  General  Equilibrium  Properties 


An  equilibrium  analysis  of  Eqs.  (l)-(3)  for  general  steady-state 

(9/9t*0)  profiles,  n*?(r),  V*?  (r ) and  <j>^(r),  proceeds  in  the  following 

J j“ 

manner.  First,  it  is  straightforward  to  show  from  Eq.  (2)  that  the 
functional  form  of  the  density  profile  n^Cr)  can  be  specified 
arbitrarily.  Moreover,  from  Eq.  (1),  equilibrium  force  balance  in  the 
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radial  direction  can  be  expressed  as  -m,V^(r)/r=e.E°(r)+e.V?  (r)Bn/c, 

3 3°  3 r 3 3°  U 

or  equivalently 


“j+ej“cj“j"ej“E“cj"°  ’ (5) 

where  e =sgne, , u)  .=  |e  |Bn/m  c is  the  cyclotron  frequency,  a>  (r)  = 

VTQ(r)/r  is  the  angular  velocity  of  a fluid  element,  and  10  (r)  is 
3 (3 

the  angular  xj^  frequency  defined  by 

“e  ” rBQ  Er  “ rBQ  3r  ^ ( ^ 

Finally,  the  equilibrium  electrostatic  potential  <J>^(r)  is  determined 
self-consistently  from  the  steady-state  Poisson  equation 

7 37  r a?  ♦ =4ne(ne-ni)  » (7) 

where  n^(r)  and  n*?(r)  are  the  equilibrium  electron  and  ion  density 
e i 

profiles.  We  note  from  Eq.  (5)  that  there  are  two  allowed  equilibrium 


values  of  uk  . Solving  Eq.  (5)  for  gives 


“j  = "j  ^ " 


■ 1 1 (' + “j  i)  ) • 


where  the  upper  sign  (w^=i0j)  corresponds  to  a "fast"  rotational 
equilibrium,  and  the  lower  sign  (10^=0)^)  corresponds  to  a "slow" 
rotational  equilibrium.  Evidently,  in  the  limit  of  equilibrium 


charge  neutrality  with  (0^=0,  It  follows  from  Eq.  (5)  that 
and  u)j=0. 


3.  ELECTROSTATIC  STABILITY  PROPERTIES 


3. A General  Eigenvalue  Equation 


To  investigate  electrostatic  stability  behavior,  we  make  use 
of  Eq.  (4)  and  linearize  Eqs.  ( 1) — ( 3)  about  the  axisymmetric  equilibrium 
state  characterized  by  n*j*(r),  ^ (^)=<*L  (r)r£0,  an<*  /3r)j£r* 


This  readily  gives 

-i  (a)- ) 6 Vj  r~  ( e j ioc ^ +2w  ^ ) 6 V..  Q 


«♦  , 

■j  31 


+tI?  (rS),rfjr  = ■ ^ “T51 

no  J 

i£rVv-o 

-i(«o-iu  >6^  + ± Ip  (rn°6Vjr)  + = 0 , 


— — r — 64> y<5<ti  = -4u£e  6n  , 

r 3r  3r  T 2r  “ j j 

r j 


(9) 

(10) 

(ID 

(12) 


where  the  index  l has  been  suppressed  on  the  perturbation  amplitudes 
5Vjr(r),  6$(r),  etc.  It  is  evident  from  Eqs.  (9)-(ll)  that  the 
perturbations  in  density  and  mean  fluid  velocities  can  be  expressed  directly 
in  terms  of  6<f>(r).  After  some  straightforward  algebra,  Eq.  (12) 
becomes  2 


I i_ 

r 3r 


i ^ 


g.  <S  4i  y 1 j) 

r j 3r 


<5  4> 


t 2 

^ <ej“cj+2“j) 

j 


2 0 2 2 

where  w , (r)=4irn.  (r)e . /m. , and  v (r)  is  defined  by 
PJ  j j J J 


Vj  (r)  = (o)-go)j ) 2-(e jO>cj+2Wj ) 


ejWcj  + r I? 


(13) 


(14) 


The  eigenvalue  equation  (13)  can  be  used  to  calculate  the  eigenfunction 


T9 


I i 


<S(|>(r)  and  complex  eigenfrequency  w for  a broad  class  of  equilibrium 


profiles  n^(r)  an^  (r)  consistent  with  Eqs.  (5)-(8). 


3.B  Eigenvalue  Equation  for  Sharp-Boundary  Equilibrium 


For  present  purposes,  we  specialize  to  the  case  where  the  electron 


and  ion  density  profiles  are  rectangular  [Fig.  2],  i.e., 


0,  \ 
n (r)  = 
e 


n =const.,  0<r<R  , 

0 p 


(15) 


, R <r<R  , 
P c 


and 


0,  , ( 0 

ni  r>  " | 


fn,,=const.,  0<r<R  , 

P 


(16) 


, R <r<R  , 
P c 


where  f=const.=fractional  charge  neutralization,  and  r=Rc  is  the  radial 


location  of  a grounded  conducting  wall.  From  Eqs.  (7),  (15)  and 


0 2 

(16),  we  find  <f>  (r)=nenQ(l-f)r  inside  the  plasma  column  (0<r<R  ). 


Making  use  of  Eq.  (6),  the  equilibrium  E *B_  frequency  can  be  expressed 

'V* 


as 


2-nen^c 

toE  = — g (1-f)  = const. 


(17) 


0 


Substituting  Eq.  (17)  into  Eq.  (8)  gives 


0)  =0) 

e e 


ce 

2 


2 (o 


yl/2) 


1 ± (l  - (1-f) 


ce 


(18) 


and 


0) 


± ci 

“i  = “i  = ‘ ~ 


/ 2“  \ 


\ I/2 


ce 


J 


(19) 


^ 2 2 

where  u>  =4Trn„e  /m  . It  is  evident  from  Eqs.  (18)  and  (19)  that  the 
pe  0 e 


angular  velocity  profiles  and  ok  are  uniform  (independent  of  r) 
for  the  electron  and  ion  density  profiles  illustrated  in  Fig.  2. 


J 


For  the  rectangular  density  profiles  prescribed  by  Eqs.  (15) 
and  (16),  the  eigenvalue  equation  (13)  can  be  expressed  as 

\ r \ 3 ' 


where 


, Mi  l 6 (r-R  ) 

r 4 2 . p 


2 2 2 
(r)=(m-£w  ) -(Ej(i)cj+2Wj)  =const. 


Vr>  ' 


(o  =const.,  0<r<R  , 

Pj  P 


0 , R <r<R  . 

P c 


2 2 a.  2 2 

In  Eqs.  (20)  and  (22),  m =47rn~e  /m  and  to  .=4irfn_e  /m.,  and  use  has 

pe  U e pi  u i 

2 a 2 

been  made  of  3<o  ./3r=-to  .6  (r-R  ). 

Pj  PI  P 

3.C  Dispersion  Relation  for  Sharp-Boundary  Equilibrium 


Except  at  the  point  r=R  , it  is  clear  from  Eqs.  (20) -(22)  that 

1 2 2- 

6$(r)  satisfies  the  vacuum  Poisson  equation  r (3 /3 r) [r36<£ /3 r ]-(&  /r  )6(j>=0. 
Therefore,  the  solution  to  Eq.  (20)  can  be  expressed  as 


6$i(r)=A(r/R  )*"  , 0<r<R  , 


inside  the  plasma  column,  and 


- / r (1"Rc  X 

6$  (r)  = A(— ) 2j> 2I~  ’ Vr-Rc  ’ 

° ' V (i-r^/r  S P 

c p 


in  the  vacuum  region  between  the  surface  of  the  plasma  column  and  the 
conducting  wall.  Note  that  the  perturbed  potential  is  continuous 
at  the  surface  of  the  plasma  column  with  6d>,(r=R  )=6d>„  (r=R  )=A. 

T i p o p 

Moreover,  the  perturbed  potential  vanishes  at  the  conducting  wall. 


U 

1 


0 
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i.e.,  <5 <(>0 ( r=Rc ) =0 . 


The  dispersion  relation  that  determines  the  complex  eigenl requency 
to  is  obtained  by  multiplying  Eq.  (20)  by  r and  integrating  from 


R (1-e)  to  R (1+e) , with  e-K),.  This  gives 
P P + 

-2  v 


T" 

3r  o 


00 


r=R 


-»  (1  - Z -I1, 


1 


3 

4r~  Sf, 

3r  i 


r=R 


(25) 


to  . e.to  ,+2oo. 

= A6$(r=R  ) 7 -f1  — 1 • 

p - v2  <o-*to 
J 3 J 


Substituting  Eqs.  (23)  and  (24)  into  Eq.  (25)  and  rearranging  terms, 
we  obtain 


0 = 


1-(R  /R  ) 
P c 


21 


l 


_Ei 


“ 2(to-£to.)  [(to-£to.)  + (e.to  ,+2<o.)] 
3 j j 3 cj  3 


(26) 


Throughout  the  remainder  of  this  article,  we  assume  that  the  electron 
and  ion  equilibria  are  rotating  in  the  "slow"  rotational  mode  with 


(0  = <0 
e e 


and 


(27) 


In  this  case,  Eq.  (26)  can  be  expressed  as 
1 


0 = 


A2 

to 

Pe 


? 0 — _ + _ 

1-(R  /R  ) 2(<o-£a>  ) [(to-lho  )-(to  -to  ) ] 

p c e e e e 


(28) 


JLL 


2 (to— «.u)±)  [(to-£toi)-(to^-toi)  ] 


+ + 

where  to~  and  toT  are  defined  in  Eqs.  (18)  and  (19).  The  linear 
e i 


dispersion  relation  (28)  can  be  used  to  investigate  electrostatic 


stability  properties  for  a broad  range  of  plasma  parameters  f, 

.2,2 

to  /to  , etc. 
pe  ce 

As  a point  of  reference,  we  note  that  the  ion  resonance  instability 


has  previously  been  studied  in  the  low-density,  low-frequency  regime 


characterized  by 


11,12 


(l-f)«l  , 


I w-iw  I <<  I Ol>  -0) 
i e i 1 e e 


Within  the  context  of  the  inequalities  in  Eq.  (29),  the  linear 


dispersion  relation  (28)  can  be  approximated  by 


0 = ± + EE 

1-(R  /R  )2* 

P c 

-2 


2(arJlcoi)  [ ( to—  Uto^)  — (oj^— u)^)  ] 


where  use  has  been  made  of  to  -ol,  and  to  -to  =co  [Eqs.  (17),  (18) 

eh  e e ce 


and  (29)].  Equation  (30)  is  identical  to  the  dispersion  relation 


obtained  by  Levy  et  al.  within  the  framework  of  a fluid  model 


in  which  the  electron  motion  is  approximated  by  V^c^xe^) /Bq . 


Throughout  the  remainder  of  this  article,  Eq.  (30)  is  referred  to 


as  the  "reference  dispersion  relation"  (RDR) . In  order  to  assess 


the  influence  of  strong  self  electric  fields  on  stability  behavior 


2 2 

[2(oj  /u)  ) (1— f ) <1  ] , it  will  be  useful,  from  time  to  time,  to 

pe  c.e 


compare  the  stability  information  obtained  from  the  unexpanded 


dispersion  relation  (28)  with  that  obtained  from  Eq.  (30).  In  this 


regard,  we  emphasize  that  the  inequalities  in  Eq.  (29)  pose  a particularly 


severe  limitation  on  the  range  of  validity  of  Eq.  (30).  The  present 


paper  is  intended  to  alleviate  this  situation  by  a thorough  investigation 


of  the  unexpanded  dispersion  relation  (28). 


P 


3.D  Dimensionless  Parameters 


It  is  evident  from  Eqs.  (17)-(19)  and  the  linear  dispersion 
relation  (28)  that  the  four  dimensionless  parameters 


m R 

f _£  —E. 

> L » „ » T> 


can  be  used  to  characterize  the  plasma  equilibrium.  Moreover,  an 
important  dimensionless  parameter  that  measures  the  strength  of 
the  equilibrium  electric  force  on  an  electron  fluid  element  is 


[Eq.  (8)] 


co_  2ai 
^ce  a) 

ce 


(1-f)  , 


which  can  be  constructed  from  the  dimensionless  parameters  in  Eq.  (31). 
Because 

2£>2 

(1-f)  < 1 (33) 

to 

ce 

is  necessarily  required  for  existence  of  the  electron  equilibrium 
[Eq.  (18)],  we  note  from  Eqs.  (32)  and  (33)  that  6 can  span  the 
range  0<6<1. 

For  future  convenient  reference,  in  Fig.  3 are  shown  plots  of 

to  /a ) [Eq.  (17)],  u± /u>  [Eq.  (18)]  and  wT/w  . [Eq.  (19)]  versus 

tj  ce  e ce  1 cl 

^ 2 2 

the  dimensionless  parameter  6 = ( 2b>  /to  )(l-f)  for  m /m. =1/1836. 

pe  ce  e i 

In  the  limit  of  maximum  charge  density  for  existence  of  the  equilibrium 

(6  = 1),  we  note  from  Fig.  3 thatto“/oo  =+(m./4m  ) and  to*  Ao  =0.5. 

I ci  1 e e ce 

Finally,  in  Fig.  4,  we  illustrate  the  region  of  the  parameter 

space  corresponding  to  existence  of  the  equilibrium 

[Eq.  (33)].  Note  that  the  equilibrium  electron  density  can  exceed  the 

* 2 2 

Brillouin  flow  limit  (to  /to  =0.5)  provided  there  is  sufficient 

pe  ce 
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2 a 2 

charge  neutralization  that  f>l-0. 5(“ce/u>pe)  . 


J 


* 

j 

\ 


4.  STABILITY  ANALYSIS 

The  growth  rate  and  real  oscillation  frequency  ior= Rew  have 

been  obtained  numerically  from  Eq.  (28)  for  a broad  range  of  plasma 

a2  2 

parameters  f,  a>  /w  and  R /R  . In  this  section,  we  summarize  the 
pe  ce  p c 

essential  features  of  these  stability  studies.  The  analysis  is 
restricted  to  nonneutral  proton-electron  plasmas  (111^/11^*1836) , 
and  the  growth  rate  and  real  frequency  are  measured  in  units  of 
the  lower  hybrid  frequency 


uLH*(<0cea>ci) 


1/2 


In  order  to  illustrate  the  need  for  an  improved  dispersion 

relation  at  large  beam  densities,  Fig.  5 shows  a plot  of  the 

normalized  growth  rate  y/wtu  and  oscillation  frequency  oi  /<oTU 

LH  r LH 

^2  2 

versus  “pg/^gg  obtained  from  Eq.  (28)  (the  unexpanded  dispersion 

relation)  and  Eq.  (30)  (the  reference  dispersion  relation),  for  £=2, 

R /R  *0.5  and  f*0.45.  Note  that  to  /toTII  is  plotted  only  for  the 
pc  r LH 

A 2 2 

range  of  “pg/^gg  corresponding  to  instability  (y>0) . Moreover,  the 

-2  , 2 

abscissa  in  Fig.  5 extends  to  “pg/^gg”©. 91,  since  physically  allowed 

equilibria  exist  for  tOpe/to^e<0.5/(l-0.45)*0.91  [Eq.  (33)  and  Fig.  4 ], 

Several  points  are  noteworthy  in  Fig.  5.  First,  the  unexpanded 

-2  2 

dispersion  relation  (28)  predicts  instability  for  0<a>pe/toce<0. 34, 

whereas  the  reference  dispersion  relation  (30)  predicts  instability 

A 2 2 

for  the  entire  range  of  allowed  equilibrium  density,  0<<Ope/ioce<0.91. 

Second,  the  maximum  growth  rate  can  be  a substantial  fraction  of  oiTU. 

Lrl 

In  particular,  for  the  parameters  assumed  in  Fig.  5,  the  maximum 


L 


growth  rate  obtained  from  Eq.  (28)  is  wt„,  which  occurs 

MAX  LH 

*2  2 

for  u /a*  "0.24.  Finally,  it  is  evident  from  Fig.  5 that  the  reference 

pe  ce  ° 

dispersion  relation  (30)  breaks  down  at  rather  modest  values  of  beam 

a2  2 

density  (oi)pe/u>ce>0. 1) , thereby  underlining  the  need  for  the  improved 
dispersion  relation  in  Eq.  (28).  Throughout  the  remainder  of  this 

section,  the  stability  analysis  is  based  entirely  on  Eq.  (28). 

*2  2 

Stability  boundaries  in  the  parameter  space  (f,  u /u  ) are 

pe  ce 

illustrated  in  Figs.  6 and  7.  In  Fig.  6,  the  solid  curves  correspond 

to  the  stability  boundaries  (y*0)  obtained  from  Eq.  (28)  for  i=l, 

and  several  values  of  R^/R^  For  a given  value  of  R^/R^,  the  region 
*2  2 

of  (f,  “pg/^g)  parameter  space  above  the  curve  corresponds  to 

instability  (y>0) , whereas  the  region  of  parameter  space  below  the 

curve  corresponds  to  stability  (ym0) . For  low  beam  densities 
„ 2 2 

(small  values  of  u /to  ),  it  is  evident  from  Fig.  6 that  the 
pe  ce 

system  is  stable  for  i,=l,  provided  Rp/Rc  is  sufficiently  small. 

On  the  other  hand,  for  high  beam  densities,  wall  stabilization 

occurs  provided  Rp/Rc  is  sufficiently  large,  i.e.,  provided  the 

conducting  wall  is  located  sufficiently  close  to  the  plasma  surface. 

In  Fig.  7,  the  solid  curves  correspond  to  the  stability  boundaries 

obtained  from  Eq.  (28)  for  R /R  =0.5  and  several  values  of  azimuthal 

pc 

a 2 2 

mode  number  l.  For  a given  density  of  physical  interest  (w^/w^^O.  01, 
say) , we  note  that  the  number  of  unstable  modes  increases  rapidly 
as  the  fractional  charge  neutralization  f is  increased  to  sufficiently 
large  values. 

The  dependence  of  stability  properties  on  fractional  charge 
neutralization  is  further  illustrated  in  Figs.  8 and  9 where  the 


normalized  growth  rate  y/ut  and  oscillation  frequency  u /w  are 

LH  IT  LH 
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*2  2 /n  2 2 

plotted  versus  f for  u>  /u>  *0.01  (Fig.  8)  and  u /u  *0.5  (Fig.  9), 
r pe  ce  ° pe  ce 

and  several  values  of  mode  number  l.  Also,  R /R  >0.5  Is  assumed  In 

P c 

Figs.  8 and  9.  In  Figs.  8(b)  and  9(b),  u /u  u Is  plotted  only  for  the 

r Lrt 

ranges  of  f corresponding  to  instability  (y>0) . Note  the  rapid  Increase 
In  the  number  of  unstable  modes  when  f Is  Increased  to  sufficiently 

large  values.  Several  important  features  are  evident  from  Figs.  8 and 

.2 

9.  First,  for  the  low-density  case  Illustrated  in  Fig.  8 (u>  /o>  >0.01), 

pe  ce 

we  note  that  only  the  fundamental  mode  (£>1)  is  unstable  when  f<0.13. 

This  Is  consistent  with  the  results  previously  obtained  by  Levy  et  al. 
a 2 2 

for  i /<i>  , f<<l.  Second,  the  maximum  growth  rate  for  each  mode 

pe  ce 

number  £ is  a slowly  increasing  function  of  f.  Evidently,  as  the 
fractional  charge  neutralization  f approaches  unity,  the  number 
of  unstable  modes  tends  to  infinity  within  the  context  of  the  present 
cold-fluid  model.  Moreover,  in  the  limit  of  equilibrium  charge 
neutrality  (f>l),  the  self  electric  field  is  equal  to  zero  and  the 
plasma  column  is  stable  (y*0  in  Figs.  8 and  9).  Finally,  for  low 
beam  densities,  we  note  that  the  real  frequency  a>r  exhibits  a 
nearly  linear  dependence  on  fractional  charge  neutralization  f 
(Fig.  8).  In  contrast,  when  the  equilibrium  self-electric  field 
is  sufficiently  large  (Fig.  9),  exhibits  a strongly  nonlinear 
dependence  on  f. 

Shown  in  Figs.  10  and  11  are  plots  of  normalized  growth  rate 

* 2 2 

y/(o  „ and  oscillation  frequency  ai  /u>  u versus  a>  /u)  for  f>0.1 
LH  IT  LH  p6  ce 

(Fig.  10)  and  f>0.8  (Fig.  11),  and  several  values  of  mode  number  1. 

Also,  R /R  >0.5  is  assumed  in  Figs.  10  and  11.  For  small  values  of 
fractional  charge  neutralization,  it  is  evident  that  the  fundamental 
mode  (£>1)  is  the  most  unstable  mode.  For  example,  for  the  f>0.1  case 


■ - - - - 
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shown  In  Fig.  10,  the  maximum  growth  rate  (ymax“0*093  ui^)  occurs  for 
o 2 

i*l  and  u /a)  *0.12.  This  Is  In  contrast  with  Fig.  11,  where  the 

pe  ce  ° 

maximum  growth  rate  for  f*0.8  (ymax=2.4  ui^)  occurs  for  £* 3 and 

,2  2 

a)  /u>  *2.5  (the  maximum  allowed  density  for  radial  confinement), 
pe  ce 

Furthermore,  for  the  f*0.1  case  shown  in  Fig.  10,  we  note  that  the 

ion  resonance  instability  is  completely  stabilized  above  some 

«2  2 

critical  value  of  beam  density  (w  /w  >0.17).  On  the  other  hand, 

pe  ce'Kj 

for  larger  values  of  fractional  charge  neutralization  (f*0.8  in  Fig.  11), 

it  is  evident  that  instability  exists  for  the  entire  allowed  range  of 
a 2 2 

u>  /w  . Moreover,  several  mode  numbers  £ are  unstable  for  a given 
pe  ce 

«2  2 

value  of  i / o»  , a feature  generally  characteristic  of  a nonneutral 
pe  ce 

plasma  column  with  significant  charge  neutralization  (see  also  Figs.  7-9), 

Of  considerable  interest  for  experimental  application  is  the 

a 2 2 

stability  behavior  for  specified  values  of  f,  oj  /w  and  R /R  . 

pe  ce  p c 

Typical  results  are  shown  in  Figs.  12  and  13  where  y/o^  and  u /o^ 

»2  2 

are  plotted  versus  mode  number  £ for  w /w  =0.01  (Fig.  12)  and 

pe  ce 

*2  2 

Wpe/^ce^O. 5 (Fig.  13),  and  several  values  of  fractional  charge 

neutralization  f.  Also,  R /R  =0.5  is  assumed  in  Figs.  12  and  13, 

P c 

and  graphical  results  are  presented  only  for  the  unstable  mode 

numbers  with  Y>0.  For  the  low-density  case  in  Fig.  12,  we  note  that 

maximum  growth  occurs  for  £=2  when  f=0.2,  and  for  £=12  when  f=0.8. 

Evidently,  for  low  beam  densities,  there  is  a broad  spectrum  of 

unstable  modes  when  the  fractional  charge  neutralization  is  sufficiently 

large.  On  the  other  hand,  for  higher  beam  densities,  fewer  unstable 

modes  are  excited.  For  example,  in  Fig.  13,  maximum  growth  occurs 

for  1*1  when  f*0.2,  and  for  i=5  when  f=0.8.  Comparing  Figs.  12  and 

13,  we  also  note  that  the  maximum  growth  rate  (measured  in  units  of  u>  ) 

Ln 

is  larger  in  the  high-density  case  (Fig.  13)  than  in  the  low-density 
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case  (Fig.  12). 

We  conclude  this  section  by  emphasizing  that  stability  properties 

also  exhibit  a sensitive  dependence  on  the  location  of  the  conducting 

wall.  This  is  illustrated  in  Fig.  14  where  the  normalized  growth 

rate  y/oj^  and  oscillation  frequency  Dr/in^  are  plotted  versus  Rp/Rc 
*2  2 

for  1-1,  (i>  / oo  -0.5,  and  several  values  of  fractional  charge 

’ pe  ce  1 ° 

neutralization  f.  Evidently,  the  1-1  mode  is  stabilized  whenever 
the  conducting  wall  is  located  sufficiently  close  to  the  plasma  surface. 
Moreover,  as  f is  increased,  the  conducting  wall  must  be  situated 
even  closer  to  the  plasma  surface  to  assure  stability.  A qualitatively 
similar  dependence  of  stability  behavior  on  Rp/Rc  Is  obtained  for 
higher  1 values. 
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5.  SUMMARY  AND  CONCLUSIONS 

In  this  paper,  we  have  examined  the  Influence  of  strong  self- 

electric  fields  on  the  ion  resonance  Instability  In  a nonneutral 

plasma  column.  The  analysis  was  carried  out  within  the  framework 

of  a macroscopic  cold-fluid  model  (Section  2) , and  electrostatic 

stability  properties  were  investigated  in  detail  for  the  case  of 

rectangular  electron  and  ion  density  profiles  (Sections  3 and  4). 

One  of  the  most  important  conclusions  of  this  study  is  that  the 

equilibrium  self-electric  field  can  have  a large  influence  on 

stability  behavior.  In  particular,  stability  properties  for  6<1 

differ  substantially  from  those  obtained  when  6<<1  (see,  for  example, 

Fig.  5).  Moreover,  for  a nonneutral  plasma  column  with  significant 

charge  neutralization,  it  is  found  that  the  fundamental  mode  (£=1) 

is  not  the  most  unstable  mode  (see,  for  example.  Figs.  8 and  9). 

However,  in  the  special  limiting  case  where  f<<l  and  6<<1,  the  J=1 

mode  can  have  the  largest  growth  rate,  which  is  consistent  with  the 

12 

result  obtained  previously  by  Levy  et  al. 

Finally,  we  emphasize  that  the  present  analysis  is  based  on  a 

macroscopic  cold-fluid  model  in  which  the  ion  (and  electron)  motion 

is  assumed  to  be  laminar.  Although  this  is  a reasonable  approximation 

when  r,.<<R  (where  r, . is  the  characteristic  thermal  ion  Larmor  radius), 
LI  p LI 

we  expect  significant  modifications  to  the  stability  behavior  when 

r,.^R  . The  influence  of  finite  ion  Larmor  radius  effects  on  the 
Li'v  p 

14 

ion  resonance  instability  is  currently  under  investigation  within 
the  framework  of  a hybrid  Vlasov-fluid  model,  which  treats  the 
electrons  as  a macroscopic,  cold  fluid  and  the  ions  in  a fully  kinetic 


(Vlasov)  manner. 
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FIGURE  CAPTIONS 


Equilibrium  configuration  and  coordinate  system. 


Rectangular  density  profiles  for  electrons  and  ions 
[Eqs.  (15)  and  (16)]. 


Plots  of  (a)  to  /to  [Eq.  (17)].  (b)  id  /to  [Eq.  (18)] 

E ce  e ce 

and  (c)  to~/to  [Eq.  (19)]  versus  6=(2<o2  /to2  )(l-f)  for 
1 ci  pe  ce 

m /m  =1836. 
i e 


a 2 2 

The  shaded  region  of  the  parameter  space  (f,to  /to  ) corres- 

2 a 2 

ponds  to  physically  allowed  equilibria  satisfying  f>l-u>ce/2u)pe. 


Plots  of  (a)  growth  rate  y and  (b)  real  frequency  to  versus 
a 2 2 

uipe^ce  obtained  from  Eq.  (28)  and  the  reference  dispersion 

relation  (30)  for  1=2,  R /R  =0.5,  f=0.45  and  m./m  =1836. 

pc’  re 

Stability  boundaries  [Eq.  (28)]  in  the  parameter  space 

(f,to2  /uj2  ) for  £=1,  m./m  =1836  and  R /R  =0.25,  0.5,  0.75. 
pe  ce  i e p c 

Stability  boundaries  [Eq.  (28)]  in  the  parameter  space 
2 2 

(f,to  / co  ) for  R /R  =0.5,  m./m  =1836  and  several  values  of  l. 
pe  ce  pc  i e 

Plots  of  (a)  growth  rate  y and  (b)  real  frequency  wr  versus 

f [Eq.  (28)]  for  w2  /to2  =0.01,  m,/m  =1836,  R /R  =0.5  and 
pe  ce  i e ’pc 

several  values  of  i,. 


Plots  of  (a)  growth  rate  y and  (b)  real  frequency  <or  versus 
2 2 

f [Eq.  (28)]  for  to  /to  =0.5,  m./m  =1836,  R /R  =0.5  and  several 
pe  ce  i e pc 


values  of  l. 


! 

f 

j 


I 

i 
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Fig.  10  Plots  of  (a)  growth  rate  y and  (b)  real  frequency  oo 

u2  /a)2  [Eq.  (28)]  for  f=0.1,  m./m  =1836,  R /R  =0.5 
pe  ce  i e p c 

several  values  of  i.. 

Fig.  11  Plots  of  (a)  growth  rate  y and  (b)  real  frequency  w 

u2  /to2  [Eq.  (28)]  for  f=0.8,  m./m  =1836,  R /R  =0.5 
pe  ce  l e p c 

several  values  of  i,. 

Fig.  12  Plots  of  (a)  growth  rate  y and  (b)  real  frequency  w 

i [Eq.  (28)]  for  w2  /u)2  =0.01,  m./m  =1836,  R /R  =0. 

pe  ce  i e p c 

several  values  of  f. 


Fig.  13  Plots  of  (a)  growth  rate  y and  (b)  real  frequency  u> 

versus  % [Eq.  (28)]  for  i2  /w2  =0.5,  m./m  =1836, 

pe  ce  l e 

R /R  =0.5  and  several  values  of  f. 

P c 

Fig.  14  Plots  of  (a)  growth  rate  y and  (b)  real  frequency  u 

versus  R /R  [Eq.  (28)]  for  <S2  /a)2  =0.5,  £=1,  m./m 
pc  pe  ce  l e 

and  several  values  of  f. 
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